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The KATRIN neutrino experiment is a next-generation tritium beta decay experiment aimed at 
measuring the mass of the electron neutrino to better than 200 meV at 90% C.L. Due to its intense 
tritium source, KATRIN can also serve as a possible target for the process of neutrino capture, 
u e + 3 H— > 3 He + +e~. The latter process, possessing no energy threshold, is sensitive to the Cosmic 
Neutrino Background (Ci/B). In this paper, we explore the potential sensitivity of the KATRIN 
experiment to the relic neutrino density. The KATRIN experiment is sensitive to a CzvB over- 
density ratio of 2.0 x 10 9 over standard concordance model predictions (at 90% C.L.), addressing 
the validity of certain speculative cosmological models. 

PACS numbers: 23.40.Bw, 14.60.Pq, 98.80.Es, 98.80.-k 



I. INTRODUCTION 

The observation of the Cosmic Microwave Background 
(CMB) stands as one of the most significant scientific dis- 
coveries of the 20th century [I, 2 . The series of measure- 
ments made after its initial observation — ever increas- 
ing in precision and scope — has helped transform ob- 
servational cosmology into an extremely predictive sci- 
ence [3]. One need only to look at the recent discov- 
eries of dark matter and dark energy to appreciate this 
progression [U |S]. Many of the parameters of the cosmo- 
logical concordance model have now been measured and 
further exploration of these parameters over the next few 
decades are likely to reveal even more hidden aspects of 
our universe. 

Though many of the predictions from cosmology have 
now been realized, there are others that remain unob- 
served. Among these is the observation of the Cosmic 
Neutrino Background (C^B) produced from the primor- 
dial Big Bang. Like their photon counterparts, their ex- 
istence in the cosmos is expected; yet, direct observation 
of CvTi remains elusive. Observation of the cosmic relic 
neutrinos — or, conversely, the absence of such — stands 
as an important verification of the theory. 

Direct observation of relic neutrinos is an extremely 
difficult challenge from an experimental perspective. 
The neutrino temperature T v is related directly to the 
CMB temperature, T 7 , such that T v = (4/11) «T 7 , or 
~1.9 K [B]. The average Cz/B neutrino kinetic energy, 
assuming a simple Fcrmi-Dirac distribution, is expected 
to be less than 0.2 meV. Most conventional methods of 
detecting neutrinos, such as water Cerenkov or scintilla- 
tor detectors, rely on interactions that have some thresh- 
old for the energy of the incoming neutrino, which is 
often many orders of magnitude larger than the C^B en- 
ergies. Fortunately, there docs exist a promising can- 
didate by which such low energy neutrinos can be de- 
tected: neutrino capture on radioactive nuclei. The pro- 
cess was originally proposed by Weinberg [7] but has 



been recently revisited by others 19] . The generic pro- 
cess, v e + N — > N' + e~ , is energetically allowed when 
m v + m,N > tojv + Tn e . Conveniently, this is also the 
required energetics at very small neutrino mass for the 
more familiar process of beta decay. The signal of this 
neutrino capture process is distinctive: the electrons cre- 
ate a monoenergetic peak at E = Eo+m^, where E is the 
kinetic energy of the detected electron, Eq is the endpoint 
energy of the beta decay process assuming zero neutrino 
mass and m v is the mass of the electron neutrino. For a 
detector with sufficient resolution, it is in principle pos- 
sible to separate this monoenergetic peak from the tail 
distribution of the beta decay spectrum. 

Recent and upcoming experiments that make use of 
beta decay nuclei have made significant advances in both 
source strength and energy resolution in order to probe 
the scale of neutrino masses. It is natural to inquire, 
therefore, what sensitivity can reasonably be expected on 
the Ci^B from these experiments given their experimental 
and physical limitations. In this paper, we explore the 
potential sensitivity of tritium beta decay experiments, 
like the KArlsruhe TRItium Neutrino (KATRIN) exper- 
iment, to cosmic relic neutrinos [TU] . 

II. DETECTION OF THE Ci^B 

The total C^B rate depends on the cross-section of 
z^-capture on tritium and the local Cz^B density. 

Rv = n„ — / o v v v f (p v ) — g- (1) 

where Na is Avogadro's number, A is the target atomic 
number, n„ is the relic neutrino density, M e g is the effec- 
tive target mass, a is the CzaB cross-section, v v and p v 
are the neutrino velocity and momenta, respectively, and 
f(pv) is the momentum distribution of the relic neutri- 
nos, which we treat as a simple Fermi-Dirac distribution 
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of characteristic temperature T v . The expected cross- 
section for relic neutrinos on tritium in the limit of zero 
momentum yields a v \ of (7.84±0.03) x 10" 45 cm 2 [5]. The 
relatively high cross-section places tritium as a preferred 
target for the sought reaction. The other commonly used 
beta decay isotope, 187 Re, is employed by the MARE ex- 
periment and is valued for its low endpoint energy. 
However, the cross-section for neutrino capture on Re 
is almost seven orders of magnitude lower than that of 
tritium. Other practical considerations, such as availabil- 
ity, the low energy of the endpoint, and final state effects 
make tritium favorable in comparison to other potential 
target materials. 

The detected neutrino capture reaction also depends 
on the local neutrino density. A local enhancement of 
the neutrino density -over the standard 55 neutrinos/cm 3 
per flavor- can occur if the neutrino mass is large enough 
to allow clumping within the Milky Way. The amount 
of neutrino over-density will in general depend on both 
the neutrino mass and the matter density profile of our 
galaxy. Table [I] shows the expected CVB capture rate 
given KATRIN's tritium target mass under a number 
of density scenarios: (a) a standard homogenous Fermi- 
Dirac distribution; (b) a Navarro, Frenk and White [T^] 
dark matter halo profile, and (c) and for a mass dis- 
tribution of the Milky Way (MW), using the local neu- 
trino densities computed in [13 . Clearly, the expected 
rates shown here are well below the detection thresh- 
old. However, understanding the analysis and limits at 
the KATRIN experiment provides a critical first step in 
understanding how to build and develop future CVB ex- 
periments. 
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TABLE I: The event rates for three different neutrino masses 
and three different mass profiles for the CzaB. Rates are cal- 
culated by scaling the results of Ref [5] by the tritium mass 
of the KATRIN experiment. All rates are given in events/yr. 



III. THE KATRIN NEUTRINO MASS 
EXPERIMENT 

The KArlsruhc TRItium Neutrino (KATRIN) exper- 
iment is the next generation tritium beta decay exper- 
iment with sub-eV sensitivity to make a direct, model 
independent measurement of the electron neutrino mass. 
The principle of the experiment is to look for a distortion 
at the high energy endpoint of the electron spectrum of 
tritium /3-decay: 

3 H ->■ 3 Hc+ + e" + v e . (2) 



The shape of the electron energy spectrum of tritium 
beta decay is determined by well-understood or measur- 
able quantities. Any deviation from this shape would 
be directly attributable to neutrino mass and would al- 
low a direct determination of the mass of the electron 
neutrino. After three years of running, KATRIN will be 
able to achieve a sensitivity of m{v) < 200 meV (at 90% 
C.L.). This level represents an order of magnitude im- 
provement on the absolute neutrino mass scale. Since 
the measurement of electrons extends beyond the end- 
point of the decay, in order to monitor background, it is 
also possible to look for the C^B signature. 

The experiment is located at the Forschungszentrum 
Karlsruhe (FZK). The FZK is a unique location as it 
hosts the Tritium Laboratory Karlsruhe, which is in 
charge of the tritium cycle of the international ITER 
fusion program. It thus offers the expertise needed to 
handle large quantities of tritium. KATRIN will use the 
windowless gaseous tritium source technique, as used by 
Los Alamos [T3] and Troitsk [T5]. Decay electrons from 
the source pass through a 10-meter long differential and 
cryogenic pumping subsection guided by superconduct- 
ing magnets. The purpose of the differential pumping 
system is to prevent gas from entering the spectrometer 
system, which would degrade resolution and raise back- 
ground by contaminating the system with tritium. 

The KATRIN experiment is based on technology de- 
veloped by the Mainz [T^] and Troitsk [T3] tritium beta 
decay experiments. These experiments used a so-called 
MAC-E-Filter (Magnetic Adiabatic Collimation com- 
bined with an Electrostatic filter) . This technology draws 
the isotropic electrons from a decay or capture event 
along magnetic field lines through a decreasing magnetic 
field so that the cyclotron motion of the electrons around 
the magnetic field lines is transformed into longitudinal 
motion along the magnetic field lines. A retarding po- 
tential is applied such that only electrons with energy 
greater than the retarding potential are transmitted to 
an electron counting detector. By varying the retarding 
potential, the shape of the decay spectrum can be recon- 
structed. The energy resolution of this measurement is 
determined by the ratio = B A — 2 oooo , where Ba 
is the magnetic field in the analyzing plane (the point 
of maximum potential) and i? m ax is the maximum mag- 
netic field. The decay electrons exiting the spectrometer 
are imaged onto a silicon PIN diode array using a 6 T 
superconducting magnet. The electrons counted in the 
detector array comprise the signal needed to reconstruct 
the beta-decay spectrum. 

The most critical aspect of the KATRIN experiment 
relevant for the detection of the Cz/B is the number of 
tritium atoms available for neutrino capture and electron 
detection. This is given by 

N efS = N(T 2 ) ■ ~ • P 

where N(T 2 ) is the number of tritium molecules, is 
the solid angle of the source as seen by the detector, 
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and P is the probability that the emitted electron exits 
of the source without undergoing an inelastic scattering 
process. 

The number of tritium molecules in the source is given 

by 

N(T 2 ) = pd-e T - A s 

where pd is the source column density, = 0.95 is the 
tritium purity, and A s = 52.65 cm 2 is the source area. 
Thus the signal rate can be written as 

n 

N eff = A s ■ e T ■ — ■ P ■ pd 
in 

The factor P ■ pd can be considered as an effective column 
density, pd e ff, which has the value 3.58 x 10 17 /cm 2 at 
KATRIN. [TO] The effective tritium source strength is 
equivalent to N e ff = 6.64 x 10 18 tritium molecules, or 
an equivalent mass of M c g = 66.5/ig. 



Because the KATRIN experiment measures the inte- 
gral beta decay spectrum above some threshold qU, the 
electron spectrum is really the convolution of the (5 and 
CVB electron spectrum, dN/dE, and the transmission 
function of the detector, T(E,qU). KATRIN also ex- 
pects a small but finite background rate, Nt, to con- 
tribute to the overall signal. Currently, this background 
rate is expected to be of order 10 mHz in the signal region 
of interest, independent of retarding voltage. 



^ —T(E,qU)dE + N b (3) 



The tritium beta decay energy spectrum has an ana- 
lytic form [T7] given by Eq. [4] 



— = Y l (NT 3 F(Z,E)p e (E + m e )(E -E)y/(E - Ef - ri*e(E -E(fs)-m v )+N c „Be ^ )P(fs) (4) 

fs 
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where F(Z, E) is the Fermi function for beta decay, E 
is the endpoint energy of the 3 H2 — > ( 3 He 3 H) + + e~ + v e 
decay, E is the kinetic energy of the emitted electron, 
and Nt 2 and Nc v b are the rates of tritium beta decay 
and neutrino capture, respectively. The gaussian term 
represents the capture signal from the Cz/B centered at 
one neutrino mass above the endpoint, with a width a 
chosen to be smaller than any characteristic resolution 
present in the experiment, but sufficiently large to be re- 
liably integrable by numerical methods. Since the target 
involves the presence of molecular T2 gas, one must in- 
clude any corrections to the endpoint energy due to the 
molecular daughter molecule following the tritium decay. 
An accounting of these states is given in 18J. Of most 
relevance are the effects of the rotational- vibrational con- 
tributions from decays to the ground state, which intro- 
duce a mean excitation energy of 1.7 eV with an inherent 
broadening of 0.36 eV. In this analysis, the final states 
are taken into account via a summation over the states 
fs of the He + T molecule, each final state weighted by 
the probability P(fs) for that state occurring. 



The transmission function, T{E 1 qU), depends on the 
value of the retarding potential, qU, as well as the intrin- 
sic resolution of the main spectrometer. For an isotropic 
source, T(E, qU) can be written analytically as: 



{0 \iE-qU<Q 
V 1 =£1 ifO<E-qU<AE 
1 if E - qU > AE. 

(5) 

where E is the electron energy, B$ is the magnetic field at 
the source, E>a is the magnetic field at qU, -B max is the 
maximum (pinch) field, and = 2 tmoo a ^ KATRIN. 
A sample decay spectrum, with and without neutrino 
capture, is shown in Figure [T] 

IV. SIGNAL EXTRACTION AND RESULTS 

To calculate the sensitivity, the data sets were fit over 
six parameters: the neutrino mass (m„), the CVB rate 
{Ncub)-, the endpoint of the beta decay spectrum (Eq), 
the beta decay rate (Nt 2 ), the width of the CVB peak 
(a), and the background rate (Nb). We assume the same 
projected KATRIN run measurement plan as reported 
in [TU]. The results of the fit in the mass and capture 
rate are plotted for 2000 simulated experiments to cre- 
ate the confidence regions shown in Fig [2] The figure 
shows only statistical errors; the systematic errors are 
estimated for four of the major errors for KATRIN: high 
voltage precisions, magnetic field precision, the effect of 
final states, and the error on the number of available tri- 
tium atoms. The errors are estimated by shifting the 
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FIG. 1: The anticipated beta decay spectrum as a function 
of retarding voltage with (black) and without (red) neutrino 
capture events. Neutrino mass is assumed to be 1 eV. 



Monte Carlo data sets by one standard deviation on the 
error and fitting at the central value of that parameter. 
The relative contributions to the errors for the mass and 
CVB rate are show in Table UU The 90% limit is shown 
for a variety of masses in Fig [3j Shown on the right hand 
side of the plot is the limit on the local density of neu- 
trinos at Earth. There is a slight decrease in sensitivity 
near 2 eV due in part because the run plan for KATRIN 
is discretized and optimized for a neutrino mass search. 
The discretization of 0.5 V in the region of interest means 
that the CVB search is restricted to a few points which 
do not change significantly with increasing mass up to 
about leV. However, above 1 eV, the finite endpoint of 
the default run plan means there are not sufficient bins 
to firmly establish the background level. Widening the 
energy scan from the original plan improves the limit sig- 
nificantly. 



Contributor 


Error (events/year) 


Percentage of Statistical 


High Voltage 


± 5850 


70.1% 


Magnetic Field 


± 2020 


24.2% 


Final States 


± 1420 


17.0% 


Normalization 


± 2080 


24.9% 


Statistical 


± 8340 




Total 


± 10680 


128% 



TABLE II: Error contributions to the CVB for four major 
KATRIN systematics at m„ = eV. Errors are extended 
to other masses as a percentage of statistical errors. Note 
that the error on the final states is limited by Monte Carlo 
statistics. 



FIG. 2: Confidence regions for cosmic neutrino captures in 
events per year versus neutrino mass in eV for four example 
neutrino masses. Statistical errors only are shown. Red ellipse 
shows 90% C.L in the Cz^B events per year and neutrino mass 
parameter space. 
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FIG. 3: The 90% confidence level sensitivity limit for relic 
neutrino over-density as a function of neutrino mass as ex- 
pected from the 3 year data run at the KATRIN neutrino mass 
experiment. Solid curves show expectation from cosmological 
prediction assuming Fcrmi-Dirac (light blue), Navarro- Frenk- 
White (violet), and Milky Way (yellow) mass distribution. 
Arrow shows neutrino mass limits already obtained from cos- 
mological observations (J] m v < 1.2 eV) [19]. 



From a limit on the local density at earth, certain the- 
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oretical possibilities can be investigated. It has been pro- 
posed that the knee of the cosmic ray spectrum could be 
caused by a relic neutrino GZK effect [5D], which requires 
a neutrino density of greater than 10 13 . Since the re- 
quired overdensity to result in the limit KATRIN can set 
is 2. Ox 10 9 , then, assuming KATRIN sees no signal at the 
expected sensitivity, KATRIN will exclude this model for 
neutrinos near earth. It has also been proposed that neu- 
trinos could couple to one another via a light scalar boson 
and form bound clouds with significant effect on small 
scale structure formation in the universe [3T] . While this 
work shows that KATRIN is able to set a better limit 
on overdensity than the experiments considered therein 
by a factor of 10 6 (resulting in an improvement on the 
limit on the fermi momentum by a factor of 100), the ul- 
timate limit on the coupling strength is also determined 
by the neutrino mass. In the range of masses accessi- 
ble at KATRIN, the limit on the coupling strength could 
either entirely rule out this model or broaden the param- 
eter space significantly. While no firm conclusions can be 
drawn with this work, it is definitely a topic for future 
analysis. 



The direct observation of the cosmic neutrino back- 
ground remains one of the strongest test of our cosmologi- 
cal framework. The KATRIN neutrino mass experiment, 
via neutrino capture on tritium, can place a limit on the 
relic neutrino density of 2.0 x 10 9 over standard model 
predictions at 90% C.L., ruling out certain speculative 
cosmological models. 
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